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OO ■ Abstract 

■ The work is devoted to the proof of the conservation of local field- 

theoretical Hamiltonian structures in Whitham's method of averaging. 
The consideration is based on the procedure of averaging of local Pois- 
on ' son bracket, proposed by B.A.Dubrovin and S.P.Novikov. Using the 

Dirac procedure of restriction of the Poisson bracket on the subman- 
ifold in the functional space, it is shown in the generic case that the 
Poisson bracket, constructed by method of Dubrovin and Novikov, sat- 
isfies the Jacobi identity. Besides that, the invariance of this bracket 
c/5 . with respect to the choice of the set of local conservation laws, used 

in this procedure, is proved. 



Introduction. 

This work is devoted to Whitham's method of averaging, which permits 
to explore the evolution of slow modulated m - phase solutions of nonlinear 
systems of equations ( see also 0, f§], @, @)- The slow modulated 
parameters of m - phase solutions (for example "running waves" if m — 1) 
satisfy in this approach to quasilinear homogeneous evolution system of type: 

Ur = Vj{U)U 3 x , i,j = l,...,N, (1) 

V=(U\...,U N ). 

For the exploring of systems of this kind it appears to be important 
to explore them on the subject of being Hamiltonian with respect to local 
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Poisson brackets of hydro dynamic type (see ||, [|J, |5J). The theory of these 
brackets, constructed by B.A.Dubrovin and S.P.Novikov, ( [|], [[|, |J), has 
been used by S.P.Tsarev (see [[7|]) for integration of systems ([I]), having 
Hamiltonian structure and reducible to the diagonal form. The investigation 
of Hamiltonian systems which do not satisfy to the last condition, was made 
by E.V.Ferapontov (see @, 0). 

In the work Q the procedure of constructing of Hamiltonian structure of 
the required type for Whitham's system of equations flU), under the condition 
that the initial system is Hamiltonian with respect to local field-theoretical 
Poisson brackets was proposed. However, the proof of the Jacobi identity 
for constructed by this way brackets was absent (see [fLOfl ). In this work 
we shall prove the Jacobi identity for constructed by Dubrovin-Novikov's 
method brackets in the case of "generic situation" with the aid of Dirac's 
procedure of restriction of Poisson brackets on the submanifold in the space 
of functions. So that, the results of this work permit to make a statement 
about the conservation of local Hamiltonian structures in Whitham's method 
of averaging. 

1. General constructions. 

Consider the evolution system on the space of fields ip = (if 1 , . . . , (p n ) of 
type: 

<Pt = Q t ( < P> ( P*> ( Px*>---)> i = l,...,n, (2) 

which is Hamiltonian with respect to local field-theoretical Poisson bracket 
of type: 

M*), ^(y)} = £ 0%{<p, <p„ . . - y) (3) 

fc>0 

(there is a finite number of terms in the sum) with the Hamiltonian: 

H = J V U (f-r, H''- (4) 

The bracket @ can be compatible with the operator of translation, that 
is, there may exist the local functional: 

P = J Vp(<p,<p x ,...)dx, (5) 
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(the momentum operator) such that: {ip l (x),P} = ip l x . (This is not a nec- 
essary condition and as was pointed to the author by O.I.Mokhov, there is a 
special class of local Poisson brackets which satisfy it, see [|TTJ, JT2]]. However, 



in the Hamiltonian structures, connected with the "physical" systems, this 
property, as a rule, is present.) Besides that, we shall admit that bracket (|3|) 
can have a finite number of annihilators, that is, functionals (not necessary 
having the form (|),(D) N U ...,N P , such that: {(p^x), N q } = . 
Definition. 

Let we are given a function of m variables $(#) = ($*(#!, . . . , m )), i = 
1, . . . , n, 27T— periodic with respect to each of the arguments, and correspond- 
ing to it m - vectors: u> = (u 1 , . . . , u m ), k = (A; 1 , . . . , k m ), such that: 

u a % a = Q i (6) 

then the corresponding to them function: 

<p(x,t) = $(kx + ut) (7) 

we shall call the m - phase quasiperiodic solution of the system (FJ). 

Note that the existence of m - phase solutions with m > 1 presents as a 
rule only for "integrable" systems like KdV. 

System (§) is a system of differential equations in partial derivatives on 
functions $(#) (9 = (8 1 , . . . ,9 m )), and its 2ir— periodic with respect to all 
variables solutions (if they exist) at all possible k and u give the full family 
of m - phase solutions of system (fj). 

We shall suppose that at any k and u in some open set (we suppose 
that system (|j) is not linear) the system @ has some family of the required 
solutions, which can be parametrized by the initial phase shifts Bq , and may 
be also by some set of additional parameters r 1 , . . . , r 9 , (the number of which 
is constant at all uj and k), which do not change under the variation of 
9q. So that, the m - phase solutions of (fj) will be given by the values 
lu = (c/,...,c/ re ),k= (A; 1 ,...,A; m ),r 1 , ...,r 9 and 9^,...,0^. 

Let now the system (FJ) has iV = 2m + g translation-invariant functionals: 

r = Jv v ((p,<p x ,...)dx, v = l,...,N, (8) 

commuting with the Hamiltonian (pj) and being in the involution with each 
other with respect to the bracket (fj): 

{r ;J P} = 0. (9) 
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Among the integrals (BI) there may be Hamiltonian (H), functional of momen- 
tum and annihilators of brackets (§) having the form (||). 

We shall assume that parameters k, u and r = (r 1 , . . . , r 9 ) are indepen- 
dent on the family of m - phase solutions of (Q) and, besides that, may be 
expressed in terms of the values of integrals (§) on the functions of family 
(0), that is k = k(U), w = w(U), r = r(U), U = (U\ . . . , U N ), where 

= & Ml - )ix = Wr jf ■ ■ f v <*' • ' 

(10) 

(it can be easily seen that values in fllCf ) do not depend upon the initial phase 
shifts 9q). So that, after the choice at each value of parameters U on M. of 
some definite function U) as having zero initials phases 9 , we can put 

values U u and 9q to be the parameters on the family of m - phase solutions 

of®. 

Each of the functionals (§) generates the Hamiltonian flow of the form: 

^ = Ql(ip,cp x ,...). (11) 

All these flows are commutative with each other and with (Q), and leave 
the family of m - phase solutions of @ invariant, generating the linear de- 
pendence of the initial phases upon the "times" r v on it and leaving constant 
all the other parameters. It means that there exist the functions cj°(U), such 
that for all $(# + 9 , U), corresponding to the family of m - phase solution 
(H), we have: 

^(U)<^=Q*(<M Q <£V,...) (12) 

From (D we can conclude that all evolution systems flllD , generated by 
integrals I u , have the same properties that the system (|2|), and all the state- 
ments, proved for system (§), are valid for systems flTTD (besides zero flows 
and operator of shift of x, we shall speak about this in the consideration of 
Whitham's method in which the nonlinearity of the system is essential). 

System (|6|) at each uj and k defines some " submanifold" (let denote it 
by -Mc^k, in the space of 2tt— periodic with respect to each 9 a functions. 
Functionals 

= ^ a % a - Q l ($, k a <$> ga , . . .) (13) 
are being the constraints, defining these submanifolds. 
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We shall assume that, if the operator I v is not the operator of momentum 
or annihilator of @, the corresponding systems (0) at all u u and k define 
(as well as (|6|)) the full family of m - phase solutions of (0), that is, the 
constraints: 

F^%\m) = - Q'M k a ^, ■ ■ •) = 0, (14) 

being consider at some definite uj v and k, define some m + g - parametric 
family ^ of 2ir— periodic with respect to each 6 a solutions, parametrized 
by initial phase shifts and some additional parameters the number of which 
is g, so that the joint of all these J^%j v k a ^ & U possible u v and k gives the 
full family of m - phase solutions of (Q) (it is so in many known examples). 

We shall also assume in our situation (it is also valid in the examples), that 
the number g of parameters r 1 , . . . , r 9 is equal to the number p of annihilators 
of (|3|), for which there exist the following motivations: 

any of m + 1 flows (|Tl|) , generated by arbitrary m + 1 integrals (|8]), are 
linearly dependent because of fll2|) on the family of m - phase solutions of 
(D, that is, there exist A X (U), . . . , A m+1 (U), such that E^ 1 ^') 2 = 1 and 
J2Y=i ^ 3 (U) = 0,Va . ^From this we can conclude that m - phase 

solutions of the system (Q) can be also defined by the relation: 

m+1 p 
j=l q=l 

(where N q - are the annihilators of (^), since the functional Y^Y^i ^(U)/^ 
generates the zero flow on the corresponding m - phase solutions), on the 
functions: 

<p(x) = $(kx), (15) 

where $(6 I ) = ^(O 1 , . . . ,6 m ) is 2ir— periodic function of 9. Supposing that 
at any k = (k 1 , . . . , k m ), fj}, . . . , /i p and A 1 , . . . , A' m+1 , satisfying the relation 
Y,™=i(X) 2 = 1, the functional J2f=i XI U] ~ Eg=i H q N q has on the functions 
(P^[) the only extremum modulo the initial phase shifts, we have that all m - 
phase solutions of @ can be characterized by 2m + p independent parame- 
ters (except the initial phase shifts), and, so that, we obtain the formulated 
statement. 

As in the finite-dimensional situation we shall require the submanifolds 
jWy^k and MqjI k (for I v which are not annihilators of @ or momentum 
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operator) to satisfy some property of regularity, analogous to the maximality 
of rank of matrix of constraints derivatives in the finite-dimensional situation. 
Namely, let us linearize the functional (0) F UM (6) = (i^ k (#), . . . , F£ M (6)) 
(respectively any of functionals (|14])) on the solution of system (]T2|), <& in {9 + 
9 , U) (it will be also a solution of flliD), that is introduce the operator Lu,6> 
(respectively L^J with the kernel L%(G + ,0' + ) (L% iu) (6 + 6 , 0' + )), 
such that: 



(2tt)» 



<^u;(u),k(u)(0) - 

1^ E L u(^ + 00, 0' + 0o)(^'(0') " $| n (0' + ^O, U))d m £ 
JO TT 



= it VA m\o) (i6) 

(similar for the constraints ^o^/kt^K^))- 

Defined by such a way Lu,e (^u^ ) are differential with respect to 6 
operators with the periodic coefficients in the space of 2ir— periodic functions 
($*(#), i = 1, . . . , n). We shall require the following conditions: 

A) For all U and 9 the kernel of operator fllq) Lu,e consists of vectors 
tangential to the submanifold A^u;(u),k(u) , defined by system (|T2|), that is, 
the functions <£v (0, k, u, r) and $ r c(0, k, u, r) give the basis in the space of 
solutions of system 

(L U)fltJ <ra)(0) = O, 

and the same property is valid for operators L^, corresponding to I v , which 
are not annihilators of (|3J) and operator of momentum. 

B) The co-dimension of the images of the operators Lu,0 o ,Lue o i n the 
space of 2tt— periodic with respect to 9 functions is equal to the dimension 
of their kernels, that is, all Lu^L^r^, have exactly m + g (g = p) left 
eigen vectors (let denote them # and «u0o )> that ^ s periodic with 
respect to each 6 a functions k^j(9 + ),j = 1, . . . , n, s = 1, . . . ,m + g and 
«w M (0 + 0o), such that: 

(^f-f *)4(M')«,0 (17) 

(2^ r- • • r ^ w (0)^ m (0, 00^0 - o. (is) 
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If all these conditions are satisfied, we say that submanifolds Aiuj,k and 
M.Ij v k in the space of 2tt— periodic functions of 9 have the property of reg- 
ularity. Besides that, under the assumption that vectors <Ev [9, k, u>, r), 
$ r c(0, k, a), r), k, u>, r) and $^(0,k, o;,r) are linearly independent at 

all values of parameters (0 o ,k, u;, r) (it is essential requirement for the fol- 
lowing consideration), the joint of described above submanifolds A4oj,k ( or 
M v UJv k) gives 3m + g = N + m— dimensional submanifold .M in the space 
of 2n - periodic with respect to each 9 a functions, corresponding to the full 
family of m - phase solutions of (0). 

2. Method of Whitham. 

The constructions described above are closely connected with Whitham's 
method of averaging for nonlinear systems of differential equations in partial 
derivatives (it can not be applied to the flows, generated by annihilators of 
@ or momentum operator), which is the following procedure: introduced 
the small parameter e, we put T = et, X = ex and rewrite the system (0) on 
the fields (p l (X,T) in the form: 

e(p t T = Q t ((p,e<p x ,(?(p xx ,...). (19) 

Let now consider the system (|19"D on the space of functions <p(9,X,T), 
9 = (9 1 ,...,9 m ) , 2tt - periodic with respect to each of variables 6 a . In 
method of Whitham we try to find the functions: 

S{X,T) = {S 1 {X,T),...,S m {X,T)), 

and 27r-periodic with respect to all 9 a functions: 

^(9,X,T,e) = ^\9,X,T,e)), 

represented by asymptotical series when e — > 0: 

oo 

¥{9,X,T,e) = E^)(^,T), (20) 

n=0 

such that the function: 
<p{0,X,T,e) = $(6+^^-,X,T,e) = fy $ (n) (#+^^, X, T, ) (21) 

6 n=0 6 
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- satisfies the system ( |I9D at all 9 and e, (e — > 0). 

It can be easily seen that the substitution of ( PH) into the system (|19|) 
gives the following equations in the zero order of e 

S$$\ 0)ea (6,X,T) =g j ($ ( o),^$ ( o)^,...), (22) 

that is, at any X and T, $( O )(0,X, T), as a function of represents one of 
the function from family .M described above, and the functions ip(x,t,e) = 
&(9 + ~S(ex, et), e), obtained from ( plj ) after the replacement of X and T 
by ex and et respectively, tend at small e to the slow modulated m - phase 
solutions of (121). 



Besides that, from ( |22| ) we obtain: 

k(U) = S x , cu(U) = S T , (23) 

where U, k and u - are parameters on A4. 
Terms with e k ,k > 0, give the relations: 

^u(x,T)Mx,T^(k)Y(e, X, T) = ^($ (0) , . . . , S x , St, - - ■) (24) 

- where Lu,e " is described in (^) linear differential (with respect to 6) 
operator with 2n— periodic with respect to 9 coefficients (expressed in terms 
of the $(o)(9, 9[o)(X, T), U(X, T)) and its derivatives with respect to 9 a ), 

is discrepancy, which depends upon the previous functions $(j)(9, X, T) and 
being of order of k, regarding that functions have order j, functions S 

- order —1, multiplication of functions adds their orders and differentiation 
with respect to T and X adds 1 to the order of function. 

The system Q2"4j ) has solutions in the class of 2tt - periodic with respect to 
9 functions if and only if at any X and T the vector ffc(#, X, T) is orthogonal 
for all defined in (|17D left eigen vectors of operator Lu(x,t),6» (x,t); that is: 

j^r C' ' ' C "SwoP + do)fl ^ x ' T)dmd = °' (25) 

which gives on the function i/,(9,X,T) at any X and T m + g = N — m 
independent relations according to the number of left eigen vectors of the 
operator L V ( X ,t),8 (x,t)- 
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Taking $ (0 ) in the form: $ (o) (0,X,T) = $ in (0 + 9q(X, T), U(X, T)), we 
can easily see from the statements formulated above that in the first order 
of e equations (2"1) have the form: 



(L[U(X,T),0 o pC,T)]$l) l (0) = 5Z a }(n)(Sx, $(0), $(O)0°'j3>(O)0«0' 3 > • • •) X 

n 

x {®U»W» u x + ®Une,e° d ox) + Fafix, $(o), $(o)*«, $(o)*«*>, • • -)Sjrx- 

- (^ (0)[ /^t + ^W^ot) , (26) 

where n = (m, . . . , n m ) is an integral m - vector with nonnegative compo- 
nents, n9 denotes here {riiO 1 . . . n m 9 m ), a!- n /3* - some definite functions, 
values Sx are St are connected with the parameters U(X), corresponding 
to $(o)(^, X), by the relations fl23|). So that the condition (p5|) gives on the 
functions U(X, T) and 6q(X, T) N — m equations of the form: 

4(U)^+BJ(U)(/1 + 4(U)^+^(U)^ = 0, ( = l,...,N-m. (27) 

Adding to them m equations: 

k° = u x , (28) 

following from (|23|), we obtain the system of N quasilinear equations on 
N + m functions U"(X, T), 9%(X, T) . 

If the conditions (p7|) and fl28|) are satisfied, the function $(i)(6*, X, T) can 
be found at any X and T from the differential with respect to 9 equation 
modulo the arbitrary linear combination of N — m functions described in (A) 
and lying in the kernel of the operator Lu(x,t)A(x,t) (let here denote them 

as £u!e )' that is ' modul ° the function of type: EsJ™ C S (X, T)^\ X)MX) . 

The values of the coefficients C S (X,T) can be found from the condition 
of the solvability of system ([24]) in the next order of e, where we have the 
same situation and, so that, if Q2"7|) and fl2"8| ) are satisfied, we can obtain 
sequentially all the terms of series (p0|), which permits to find the required 
functions $(6*, X, T, e). Function S(X, T) can be constructed using the values 
k(U) and w(U). 

Lemma 1. 
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Under the conditions formulated above, that is, in the presence of N local 
translation- invariant integrals @, such that the evolution of their densities 
V v (ip, <p x , . . .) according to the flow (0) has the form: 

7? = K£(¥>, ¥>„...) (29) 



with some 1Z U , the system ( P7|) does not contain o (X,T) and provides the 
relations only for U(X, T), that is, all the terms containing 9qx and 9qt are 
orthogonal to the left eigen vectors of the operator L independently upon 
the values of parameters XJ(X,T) and 9 (X,T), so that A'j* = 0, B'f; = 0. 
Besides that, in method of Whitham for any of the systems (|TT|), which is 
not the operator of translation or trivial flow, we have the same situation. 
Proof. 

If (f?7]) and ( JZSP are valid, there exist the solutions of (|T^) in form of 
asymptotic series (|20|), the substitution of which into fl29|) and integration 
with respect to 9 gives in the first order of e: 

UZ = d x {K v )(U), (30) 

- where (...) means the averaging on the family J\A defined by the formula: 

I i-2n r2n 



(f(<p, <p x , . . .))(u) = — — r. . . r f($, k a ^,k a k^ eaeB , . . .)d m e, 

(2n) m Jo Jo 

(31) 

if (p(x) = <5> in (kx + 9 ,U). 

So that, from the system follows the system (p0|), which is a 



system of independent in the general case equations on functions U(X, T), 
having the same form as (p^)-([28|), and, so that, (|30| ) is equivalent to (p7|)- 



!8[). Since these conclusions may be applied to any of systems (0), which 
is not translation or trivial flow, the Lemma is proved. 

Let us note that in the presence of the additional conservation laws of 
form fl2"§D their averaging gives the equations which are the corollaries of 



(|30D, the independence of (23), which is the system of N equations on iV 
parameters U U (X), is the property of generic situation. 

System ( |30D (or equivalent to it (p7|)-(p8D), giving the evolution of slow 
modulated parameters U(X) of m - phase solutions of (§) is called the system 
of equations of Whitham. 

System (|30|) has the form: 

UZ = VZ(\J)U£, v,n = l,...,N, (32) 
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that is, it refers to evolution systems of hydro dynamic type. For smooth 
initial dates U(X) system (|32| ) has in general case smooth solution up to some 
moment T depending upon the initial dates and after that the solution will 
be broken. So that we can use the solutions of Whitham up to the definite 
point of time, after which they are not defined. 

Let us point also that usually we consider just the finite number of terms 



of asymptotic series (|20|) when they give us solutions of initial system (0) 
modulo terms of higher order of e. In particular it is possible to consider just 
the first term of (|20|) , which gives the evolution of slow- modulated m-phase 
solutions of (0) provided that the Whitham equations on parameters U hold 
and the next term in (^) is uniformly bounded. We shall consider all arising 
here asymptotic series from this point of view, in particular, we shall not 
interest in there convergency regions since we consider the relations on the 
first (say k) terms of such series modulo the higher orders of e|] 

3. The Poisson brackets of the hydrodynamic type. 



Among the systems (p2|) there is an especial class of them which are 
Hamiltonian with respect to Poisson brackets of type: 

{U U (X), U»(Y)} = g ull (U(X))6'(X -Y) + b^(V(X))U^5(X - Y) (33) 

with local Hamiltonian of hydrodynamic type: 

H = J h(V(X))dx, (34) 

which plays an important role in their integrability (see |7|], 0, 0). 

The theory of brackets (|33|) with nondegenerated g uti (XJ), constructed by 
B.A.Dubrovin and S.P. Novikov (see 0, 0, 0), is closely connected with 
Riemannian geometry. In particular, from their skew-symmetry follows: 

<r = <r, bT+K v = ^, (35) 

Leibnitz identity leads to the fact that under the transformations of coordi- 
nates U v — > U U (XJ) functions g vtl transform as the contravariant components 



1 Author is grateful to I.M.Krichever for fruitful discussions of these questions. 
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of a metric tensor, whereas the functions = —g^ T b\ u (g^ T g TU = 5") trans- 
form as coefficients of connection, consistent with metric view (051) • The 



Jacobi identity for ( ]3"3| ) in the case of nondegenerated metric g u ^ is equiv- 
alent to the symmetry of connection and zero curvature of the metric: 

The theory of brackets (|33|) with degenerated g Ufl (XJ) is more complicated, 
but also has a nice geometrical form, see |Tj| . 

In Q B.A.Dubrovin and S.P.Novikov proposed also the method of con- 
structing of the brackets fl33]) for Whitham's system of equations (H), start- 
ing from the Hamiltonian structure (Q) for the initial system (|2|) under the 
condition that we have the necessary number of commuting integrals (§). 
Namely, let us calculate the brackets of densities of integrals (|§p in the form: 

{V(<p(x), <p a , . . .), P"(v>(v), ¥>„...)} = £ A?{<p(x), <p x , . . .)5\x - y), (36) 

fc>0 

(there is a finite number of terms in the sum, v, /j, — 1, . . . , N). 
View (H) we can conclude that 

A^{y,y x ,...) = d x Q^{y,y x ,...) (37) 

for some Q vtl . In the described above coordinates U = (U 1 , . . . , U N ) the 
brackets of Dubrovin-Novikov have the form: 

{U»(X),U>*{Y)} = {AT) (U{X))6'{X - Y) + 9 <^>(J J W) S ( X - Y) (38) 

oX 

-where (...), as previously, denotes the averaging on the m - phase solutions 
of (0), defined by formula (|3~ID . 

However, the proof of the fact that constructed by such a way brackets 
satisfy the Jacobi identity, was absent in [|J (see |H| ). The main purpose 
of this paper is to prove the fact that the procedure (|36| ) - ( p8|) of averaging 
of brackets (|3|) really gives the Poisson brackets of type (|3"3"D, satisfying the 
Jacobi identity. 

We shall need for the further purposes the Dirac procedure of restriction 
of Poisson bracket on the submanifold. Let us describe here this procedure 
in the notations of spaces of finite dimension. 

Let in the space V with the coordinates x = (x 1 , . . . , x ) and the Poisson 
bracket: 

{x 9 ,x p } = J 9P (x), (39) 
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be a submanifold N, defined with the aid of constraints g , . . . } g s by the 
equations: 

ff 1 (x) = 0, ... ,/(x) = 0. (40) 

Let us suppose also that some I — s functions / 1 (x), . . . , /'~ s (x), defined in 
the vicinity of N in the space V, give the coordinate system on N after the 
restriction on it. All the redefinitions of the functions / A (x), having the form 

/\x) = /\x) + £r c W(x) 

with arbitrary t a (x), do not change this coordinate system. Let us suppose 
that we can find functions t^(x), such that the submanifold N is invariant 
under Hamiltonian flows generated by all the functions / A (x) (in the case of 
nondegeneracy of the matrix {g^(x), <r(x)} on N it is always possible), that 
is: {/ A (x), g^(x)} = if g(x) = 0. Then we can define the Dirac restriction 
of bracket ( 09|) on the submanifold ( |4~0|) by the formula: 

{/ A ,n* = {/ A (x),/>(x)}u(f) (4i) 

The bracket (O) satisfies automatically all the necessary identities. It can 
be easily verified that with the functions r A (x), founded by such a way, the 
bracket fl4"T| ) may be also written in the form: 

{f\rr = {/ A (x),r(x)}u(f) - (^(x)^(x){^(x)^(x)}) \ N (f). m 

The infinite dimensional form of (^) will be very convenient in the further 
analysis. 

4. The coordinates in the vicinity of submanifold corresponding 
to the full family of m-phase solutions. 

Let us now return to the constructions connected with Whitham's method 
of averaging. After the substitution X = ex the bracket @ has the form: 

{<p\X), ^(Y)} = J2 B%(cp, e<p xx , e W • • ■)e k 5 ( - k \X - Y), (43) 

fc>0 

integrals (||) become: 

F = e- X Jv(<p,e<p x ,...)dX (44) 
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Let us consider the space of functions ^(O 1 , . . . ,9 m ,X), 2n— periodic with 
respect to each of 9 a and define at all e the Poisson bracket on it according 
to the formula: 

{<p% X)^{9\ Y)} = £ Bi j (<p, e Vx , eVxx, • • -)^ k \X - Y)5(9 - 9') 

k>0 

(45) 

(where 5(9 — 9') is the 5— function in m— dimensional space). 

Consider in the space of 27r— periodic with respect to 9 functions the 
submanifold Ai' of functions if (9, X), such that at any X ip(9, X) as a func- 
tion of 9 lies in Ai. The functions U U (X) and 9q(X) can be taken as the 
coordinates on the submanifold Ai', so that the functions (p(9, X) from Ai' 
will be represented by the formula: 

<p{6,X) = <5> in (9 + 9 (X),U(X)). (46) 

After the prolongation of coordinates U(X) and 9q(X) by the independent 
upon e way in the vicinity of Ai' (let us denote it by A$) in the functional 
space of 27T— periodic with respect to 9 functions, where these functions 
satisfy the conditions: 

(p(6,X) e A 8 & 3 Vo (9,X) e M' : \\<p{9,X) - ip Q {9,X)\\ = 
= ( max I <P* (9, X) - ifl (9, X) I + — (J2 max \ ip l 0a - ifl 0a \ + max | tp x - <p % ox | ) + 

i ^- a 

+ ^ (J2 max I Ve-e/s - ¥oe«,en I + J2 max I ~ ¥o0<*x I + max\(p xx - ipl xx | ) + 

(47) 

+!(...) + ...)<« 

(that is we imply that there exists a function <p Q (6, X) from Ai' such that all 
the terms of the series (f4"T|) exist, the series fl4"?|) is convergent and satisfies 
to the formulated condition), we can define the submanifold Ai' with the aid 
of constraints: 

F*(9,X)m = a; Q (U[^(l))^ - Q\ Vl A; Q (I%](X))^, ■ ■ ■) = 0. (48) 

The difference between the systems (U) and fl48|) is that the system does 
not depend upon the parameters k and u, since they are now determined 
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functionals of ip(9, X) (in the vicinity of M.'), the values of which on the Ai' 
coincide with the corresponding parameters of the family. It can be easily 
seen that only the functions from J\4' satisfy to ([f5|). 

The functionals U(X) and 9q(X) i n can be defined for instance by 
the following way: let introduce in As N functionals of the form 

a>](X) = -i- r... r A v { V {d,X),y ea {e,X),...)<re (49) 



(27r) m Jo "Jo 

with some functions A v , such that on the M! the values a u are functionally 
independent. On the .M' the values a v , as can be easily seen, can be expressed 
in terms of the values U(X) and do not depend upon 9q(X). So that, dividing 
M' into the "maps" in which \UKJX) - U? 2 JX)\ < 5', we can express U(X) 
in terms of a u (X) in each of these maps in the form U U (X) = f u {a.{X)) and 
then, using the definition (f49|) of a(X) in As, we can extend U"(X) into the 
vicinity of each map (this can be done at each X independently). Then, 
after the definition of the functionals U U (X) in As, let consider in As the 
functionals: 

I p2tt /*2"7r 

^W = ^y o ■■■J o E^(^)C(e,%]W)«n, (so) 



where U) are introduced in ( ^6l) functions from A4'. 



If <p(0,X) G M' and (X) = 0, then $ a (X) = and, in the generic 
situation, at small 0q(X) the values {9q(X)} can be expressed in terms of 
{iT(X)} on M! in the form: 6g(X) = T a (i3(X)). After that, by the same 
way, dividing, if it is necessary, each of the maps described above into the 
parts in which: |0o(i)PO — ^o(2)POI < (independently at each X), and 
expressing by such formulas 0q(X) in terms of the $(X), we can extend 
Qq{X) into As using the functionals $(X). 

It can be easily checked that after such definition of U(X) and 9q(X) in 
As we have for two functions ip^(9,X) and ip, 2 -\(9,X) from As, satisfying 
the condition of type (^7[), that is \\(p^(9,X) — ipr 2 \(6,X)\\ < 5, where 5 is 
small, the relation: \Ufa(X) - Ufa{X)\ < C5,\9% (1) {X) - 6% {2) (X)\ < CS, 
where C is constant. The analogous relations will be also satisfied for the 
variational derivatives of U(X) and 9q(X). 

We shall suppose that after the prolongation described above the sub- 
manifold M! ', given by the system (f4"gp, possesses the property of regularity, 
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analogous to the regularity of M.uj,\l- Namely, it is clear that the vectors: 
§ ea {0 + 9 {X),V{X)) and <5>uu(e + 6 {X),U{X)) (let us denote this set as 

{£iu 0i (X) , q = 1, . . . , N + m}) lie at all X at the kernel of linearized on 
the function from Ai' (characterized by the corresponding values U(X) and 
9 (X)) functional F i (9, X)[tp], introduced in ([48]): 

5F\e,X) = ^^^ 2 ^.^L} [XJ] (^o(X),^ (X),X)V(^^)^' = 



2vr i-2-k 



r zn dO i 
J [(u a (U(X))5}8 e «(9 - 9') - -^-\ k=const 5(9 - 9') 



{2n) m Jo Jo LV v v " 3 y ' dip- 
8Q i dQ i 

— — j \k=amst&6 a \9 ~ 9) — ——■ |k=consi^0«0/3 {9 — 9) — . . .) + 

du° 6U»(X) d® 6U»(X) 

+ { d!F {x) s^(9',xf ea ^ x) ~ oW^ x) WW^) )] * ( ' } ' 

(51) 

We shall assume that there is no other vectors possessing this property 
and, besides that, at all points of M! there exist at all X exactly iV + m of 

"left eigen vectors" /{ ^ i(Jf) for the operator L[u,e ] with zero eigen values, 
that is such linearly independent at each X 2ir— periodic functions of 9: 
N ( [i](# + #o,X), i = l,...,n), that: 

1 f\ . . fj ~k% (9, X)U m (9, 9', X)d m 9 = 0. (52) 



(2?r)" 

Note that unlike (k%,e }{X)} the form of functions {k\^j e ^(X)} depends 
upon the manner of prolongation of U(X) and 9 (X) into the As, since 

from the kernel of L correspond to the variations of functions (p(9,X) 
along the J\4', on which U(X) and 9q(X) are defined initially, whereas the 

fS q \9,X) are connected with the image of the operator L, for finding of 
which it is necessary to know the change of U(X) and 9 (X) at all variations 
of <p(9,X). 

With the aid of the constructions described above we can introduce an- 
other system of coordinates in As instead of the standard system, given by 
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the functionals ip l (0, X). Namely, we take as the coordinates in As the func- 
tional {U V {X) : 9%(X), v = 1, . . . , N, a = 1, . . . , m}, and 

G\ UM [ip] (9, X) = — L; fj. . . V m (9 + 9 (X), 9' + 9 (X), X) x 

x {^{9',X) -<% n (9' + 9o(X),V(X))d m 9', (53) 

a i 

where Ljnji are functionals from ([51]), appeared in the linearization of the 
constraints F l (9,X) on the submanifold Ai'. 

So that, the functions ip(9, X) from will be characterized by the set of 
N + m smooth functions U"(X), 9q(X), and by the functions G l (9,X),i = 
1, . . . ,n, taking the values at the space of 2ir— periodic with respect to 9 
functions, satisfying at each X (and given U(X) and 9q(X)) to N + m linear 
integral conditions of type: 

fj. . . fj k% (9 + 6 Q (X), X)G i (9, X)d m 9 = (54) 



(27T) 

(The conditions of such type are not customary in usual systems of coordi- 
nates, however, they can be encountered in the theory of stratifications). 
Lemma 2. 

At small enough 5 the values of functionals U u (X),9q(X) and G l (9,X) in 
A§, satisfying (|5l]), give uniquely the values <p l (9,X). 

Proof. 
Indeed, the system: 

(2^ c ■ ■ r ^ + uxi °' + ux)ix) x 

x{^{9',X) - + 9 (X),XJ(X))d m 9' m = G\9,X) 



has, view (p4|) , a solution, defined modulo the linear combination of the 
vectors {£[u0 o ]PO}j tangential to M.' and corresponding to variations of 
parameters U(X) and 9q (X); so that, their coefficients in small enough A^ 
may be defined by the values U U (X) and 9q(X). 

We shall also need another system of coordinates in As, which connected 
with the system described above by the transformation depending upon e. 
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Namely, consider the functionals: 



J V {X) = j^- J*. • • f* V&(9, X), ep x (9, X),.. .)d m 9, v = 1, . . . , N, 

(55) 

and 

e*(x) = e (x) - - [ x k(j(x'))dx' (56) 

(for some initial point X ). As a system of coordinates in the vicinity of 
M! we take the functionals {J u (X),9q(X)}, and also the set of constraints 
&(9,X) satisfying to (0). 

The transition from (U U (X),6$(X)) to ( J V {X), 9* a (X)), as can be easily 
seen from ([SB]), is defined at e 7^ 0. On the submanifold JW' (that is at 
G*(0, X) = 0) at e -> we have: 

^ W = ^ f . . . /V(* ta (* + W, u ( x)), 

= e J(fe)(U, Ux, • • • , Ujtx, $ox, • • • > Qokx), 

fc>0 

where JY k s represents the averaged with respect to 9 the fcth Teilour term 
of the expansion of V v at e — > on the functions from .M'. There is a 
finite number of terms in the sum, there introduced the notations: \Jkx = 
d fc U/<9X fc ,# fcx = d k 9 /dX k . The expression for 9*(X) in terms of U(X) 
and 6*o (X) on the J\4' can be obtained by the substitution of these expansions 
for J(X) into (|56|). 

We shall need also the inverse transformation from J(X) and 9q(X) to 
U(X), 6> (X) at G q (9, X) = (that is on the M'). Let us note that the values 
J^(X),6?Q a and U^(X) are connected on M.' by the relations (the definition 
of J V {X)): 

r(X) = j-L- \ 2 \.. rV(^ n (9 + 9* (X) + - / X k(J(X'))rfX',U(X)), 
(ztc)" 1 Jo Jo e Jx 

ed x $ in {9 + 9*{X) + - [ X k(3(X'))dX', U(X)), . . .)d m 9 = 

€ JXo 
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= 7^w f...fr(Mff+^ + ) / x k(j(x'))rfx', u(x)), 

{2Tc) m Jo Jo e Jx 

k a (J)d e ^ in (9 + 9*(X) + - ( X k(3(X'))dX', U(X)), . . .)d m 9+ 

e Jx 

% e *WrL -L "&)<•*(■■■>,..)«, 

- where ••)>•••) are l° ca l functionals depending upon Qi n (9+9 Q f (X) + 

± f£ k(J(X'))dX',U(X)) and their derivatives with respect to U" and 9 a 
with the coefficients of type: XJx(X), TJxx(X), ■ ■ ■ , k(J), 9xk(J), 
<9|-k(J), . . ., and 6>q X (X), . . ., given by the collecting together these 

terms, having the general multiplier e k . The term, corresponding to the zero 
power of e, is written separately. 

After the integration with respect to 6, which removes the singular at 
e — > phase shift 9q in the argument $j n , we obtain on J\4': 

J V {X) = C(J,\J)+Y^ e fc C(fc)(U, Ux, • • • , Ufcx, J, Jx, • • • , Jfcx, 9* 0Xl . . . , #ofcx)- 
fc>i 

(57) 

The sum in (|57|) contains the finite number of terms, functions (fy and C 
are integrated with respect to 9 functions VY k \ and V v respectively. Since 

u) = r ... r v v {® in {9, u(x)), ^(j)^$ m ^, upo), 

(Z7TJ m JO JO 

fc a (J)^(J)^ Q ^$ m (^, U(X)), . . .)<T0 

(in the integration with respect to 9 we omitted the unessential phase shift 
$, n ) we obtain that the system: 

J^(X) = C(J(X),U(X)) (58) 

is satisfied by the solution J V (X) = U U (X), according to the definition (|T0|) 
of parameters U(X) on the family M.' . Since we suppose that the system 
fl5"8"|) is of general form, we shall assume that it is equivalent to the system 
J V {X) = U V {X). 

Taking this into account, we can resolve the system fl57p by the iterations, 
taking on the initial step U V (X) = J U (X)\ for U V (X) we shall obtain the 
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expression of the form: 



u»{x) = r{x) + £ e k u( k) (j, j x ,..., j kX , e* x , e* kX ). (59) 



The substitution of ( |59"D into (]57f) , under the condition of nonsingularity of 
matrix || ^J§^-\\ |u=J> wm uniquely define the functions UY k y The value 
9q(X) is expressed in terms of J(X) and 9q(X) according to the formula 

Later it will be convenient to write the expressions of type 9q[J, 9q, e] and 
U u [J, 9q, e], assuming the expressions (|56| ) and (|9|) (not necessary on the 
M). 

The functionals (|55D are well defined on the full space of functions f(9, X) 
and the Hamiltonian flows, generated with the aid of (|45|) by the integrals: 



J q (x)r(x)dx = J q(x)-^— j^...J^v v ( v ,e^ x ,...)d m 9, (eo) 



(27 

have for all smooth functions q{X) the form: 



$ = Q(X)Ql(<P, ei Px, ■ ■ •) + £QxQl(<P, ecp x , . . .) + e gxxQ^, e^x. •••) + •••> 

(61) 

where QJ, are introduced in (pTTj) flows generated by the functionals I u , Q l u , 

x i 

Q v , . . . are some functions. (The derivatives of q(X) with respect to X arise 
in the calculation of the variational derivative and also as a result of action of 
bracket (^). As can be easily seen, each differentiation of q{X) with respect 
to X appears with the multiplier e.) 

5. The restriction of the Poisson bracket on the submanifold 
corresponding to the full family of m-phase solutions. 



For the Dirac restriction of bracket P5| ) on the submanifold M.' in the co- 
ordinates 3(X),9* (X) and G(9,X) = (G l (9,X)) we shall need their Poisson 
brackets with each other (on the M.'). 

We shall begin with the brackets of type: {J V (X), J»{Y)}. View (g|), 
these brackets have the form: 

{r(x),j^Y)} = 
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l Z7r J J0 fe>0 

/*... / a, E^%,-)^ ) (^-w (62) 

where view (P7|): 

^o M (^, eVxx, • • •) = ed x Q vtM (<p, eip x , e 2 ip xx , ...). (63) 

At the points of the submanifold M.' the function (p(9,X) has the form: 
V \9, X) = $1(9 + 9*(X) + i /* k(J(X'))dX', U[J, 0J, e](X)) and after the 
substitution of it into ( |6"2|) we can obtain the brackets {J U (X), J^{Y)} at the 
points of M.' (characterized by coordinates J(X), 9q(X)) in form of regular at 
e — > series, since the integration with respect to 9 removes the irregularity 
at e — > in the argument of function $j n . View (|63"D, it can be easily seen 
that the zero terms are absent in these series, so that they can be written in 
the form: 

{r(x),j»(Y)}\ G(f):X)=0 = 

= e ((A7)(J(X))5'(X — Y) + d{QV ^ X)) 5{X - K)) + e 2 J^[3, 91 e], 

(64) 

where J ViX are regular at e — > functionals of J and #q. The functionals of 
this type will arise later and have the form of regular at e — > asymptotic 
series of type: 



k 



C[J> ^o> e ] - C(fe)(J> Jx, • • • ) JfcX, 6*0X5 • • • J ^0fcx) e 
fc>0 

Here we used the relation ( |59|) giving U(X) in terms of J(X) and #oPO 
at the points of A4', according to which the arguments U U (X), on which 
depend the values (A^) and (Q"^) on Ji4', are replaced modulo the terms of 
higher order of e by J U (X), (. . .) means, as previously, the averaging on the 
family of m-phase solutions of (0), defined by the formula (|3lQ . 

The evolution of the densities of conservation laws V u ((p, €ip x , . . .) accord- 
ing to the system (|6l"D, generated by the functional / q(X) J^(X)dX, as can 
be easily seen from (|36|) , has the form: 

V^X) =^A?(( P ,w xx ,£ l <p xx ,...yq kx (X). (65) 

fc>0 
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The evolution of J U (X) at the points of A4' , characterized by the coordinates 
J(X), 8q(X), will be written in the form: 

J V AX) = e ((A7)(J(X)) qx (X) + 9( ^ J(I)) g(I)) + 

+ e 2 f^[J,9*,q{X),e), u, fi = 1, . . . , N, (66) 

where f Ufl are regular at e — * asymptotic series with respect to e k , the 
coefficients of which are local functionals of q(X), 3(X) and 6q(X). 
Lemma 3. 

From the system (|66|) [y = 1, . . . , N, \i - is fixed), giving the evolution of 
J(X) according to the system (|6T|), at the points of M.' follow the relations: 

k%(J (X)) = e (q(X)u«(J(X))) x + e 2 ^[J, 6* , q(X), e], (67) 

- where lu u = (ul, . . . , uj™) are introduced in (0) functions, are regular 
at e — > functionals of J, 9q and q(X) (linear with respect to q(X)). 
Proof. 

System (|6l|) admits the construction of the family of solutions similar to the 
described in Whitham's method ones, that is the family of solutions of type: 

<P = V {^f 1 + 0, X, T, e ) = £ Vk (^fl + 9, X, r) e\ (68) 
T = er^. 

By the substitution of (^) into ([H]) we obtain in the zero order that 
ip (6,X, T) at each T is the function from the Ai f , since it satisfies to one 
of the systems (|l"2"D (we assume here that J M is not the momentum operator 
or annihilator of bracket (^), for which the statement of Lemma is evident), 
besides that: 

of = q(X)u«(\J(X)), a\ = k a (XJ(X)). (69) 

We assume, as it was declared previously, that the systems (|TT|), generated 
by the funcionals I** (except the annihilators and momentum operator), and 
corresponding to them submanifolds -Mo; Mi k satisfy to the nondegeneracy 
conditions (A) and (B), which permit to construct the asymptotic series 
fl68|). The conditions of the solvability of system on (fi(9, X, T), that is N — m 
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equations of orthogonality of the discrepancy fx to the left eigen vectors of 
the corresponding operator L[u(X),0 o pO]> together with the equations 

k(U) T = (q{X)ui»(U)) x , (70) 

following from (|69"D , give, as previously, the system of equations on U(X) 
and 6 (X). 

It can be easily seen (view (|69|)), that the system on $(i)(6*,X, T) has 
now the form: 

q(X) (tfa {x)MX)] $x) (0, X, T) = g(X)E afflux, *<D), *(0)«-, ■ • 0* 

n 

X ($(O)n0,t/" 

+g x (X)Q;($ (0) ,a x $ (0) , Q , . . .) - (Sfo^O? + ^flSr) , 

all the notations are the same that the notations in (|26|), the values <tt and 
a x are connected with the parameters U(X), corresponding to <&(o)(#, X, T"), 
by the relations (|69|) , the functions (XjA and are the same that in 

Whitham's method for the evolution system, generated by the operator I v . 
The last property permits to conclude that, similar to Whitham's situation, 
the conditions of orthogonality of discrepancy fx to the left eigen vectors of 
L do not impose any restriction on the parameters 9q(X, T) and have a form: 

q(X)Ai(V)U x + Bi(\J)U» + q x (X)C<(U) = 0, ( = l,...,N-m. (71) 

Together with the equations (|70|) , the system (|7l|) gives N equations on the 
functions U(X). 

At given initial dates the system (f70|) -( [7T|) has the smooth solution at 
T < T (up to the moment of destruction of waves), and in this region, as 
in Whitham's method, we can, using the function (p, Q J9, X,T) from Ai', 
characterized by the values U(X,T) and 9 (X,T), construct uniquely the 
next terms ipr k \(9,X,T) of the series (|68l), which will be the functionals of 
U(X,T) and 9q(X,T). Taking this into account, and using the definition of 
J(X) ( p^D and the parameters U(X) ([H]), we can obtain on the family flB5|) 
of solutions of (^1|) the relations: 

r(X,T) = U V (X,T) + ^2e k J^ k) [U,9 ], (72) 

fc>i 
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- where U(X, T) and 6q(X,T) are parameters, connected with the function 
<yj( ) (9, X, T) from M.', J(X) are values of the corresponding functionals on 
the solutions (|B8D, corresponding to the given XJ(X,T) and 6*o(X, T). 

After the substitution of ip(9,X,T) in form of ( |S8|) into the expression 
( 65]) and the integration with respect to we obtain the expression for the 
evolution of J(X) on the solutions fl68|), which, as can be easily seen, has the 
form: 

j; M (X) = 6(g x (X)(^)(U(X)) + g(X)9 x (Q^)(U(X))) + 6 2 r At [U,0 o ,e]. 

(73) 

Comparing flTSQ and ( |72|) we can conclude that the system 

U^ = q x (X)(Ar)(V(X)) + q(X)d x (Q^)(U(X)), (T = er M ), - 

is equivalent to (|70|j7l|) , and so that, the relations (|67|) follow from (|66|) . The 
Lemma is proved. 
Corollary. 

On the functions of family A4', characterized by coordinates J(X),8q(X) at 
e — > 0, take place the following relations: 

{A;«(X), J^(y)} = e ^(J(X))5'(X - F) + ^^ X)) 6(X - V)) + 

+ e 2 £^[JA*,e],- (74) 

where /c a/i are regular at e —>■ functionals of J(X) and 8q(X). 
Lemma 4. 

Under the formulated previously assumption that the number of parameters 
U u is equal to 2m+p, where p is the number of the annihilators of the bracket 
@, the following conditions take place: 

E^P^(U)-0, Va,/?. (75) 

!/=l U 

Proof. 

As was formulated previously, any of m+1 flows, generated by m+1 integrals 
from the set @, are linearly dependent on the family of m - phase solutions 
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of system (g), and for some Ag(U), v = 1, . . . , N, s — 1, . . . , iV — m we have 
on Ai' (N — m) independent relations: 

N 

E^(UK(U) = 0, a = l,...,m, (76) 

u=l 

- where u;"(U) are introduced in (|12"D frequencies, corresponding to the inte- 
grals I u . It follows from this that on m - phase solutions of (§), characterized 
by the parameters XJ(X), for some functions /^(U) take place the relations: 

N P 

Y,K(V)6I[cp] = Y,» q s(V,0o)W q , - (77) 

u=l q=l 

where N q are annihilators of bracket @. Let us denote the values of the 
annihilators N q on m - phase solutions of (^) by n 9 (U) (the dependence of 
#o is absent because of the commutation of N q with I", generating the linear 
dependence of the initial phases upon the time). The values of the functionals 
Fonm - phase solutions of (0) are the parameters U u , and for infinitely small 
shift on the M. along the direction £ (in the space of parameters on the Ai), 
tangential to the submanifold k(U) = const, (k = (k 1 , . . . , k m )), we obtain 
the uniformly bounded variations of functions (f l {x), so that, view ([77]), we 
obtain on M.: 

N p 

£ XZMdU" - £ a*2(U, Oo)dn q (V) = 0, 

v=\ q=l 

if dk a = 0. 

From this we can conclude that /i(U, 6*o) = A*(U) (the dependence upon 
6*o is absent), and for some functions r"(U), a = 1, . . . , m, s = 1, . . . , iV — m, 
take place the relations: 

N p m 

Y: KWdU" = £^(U)rfn q (U) + ]T r s «(U)^(U). (78) 

v=l q=l a=l 

After the expression of iV — m = m + p differentials dn q (XJ) and dk a (XJ) from 
N — m equations (|78|) (we assume that in the generic situation this can be 
done) we obtain the relations: 

N-m N 

dk a (V) = £ a as (V) £ A^(U)d^, (79) 

s=l y=l 
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N—m N 



dn q (XJ) = &S(u) E K(V)dlT (80) 



8=1 V=l 



for some functions a as (XJ) and &*(U). After that the statement of the Lemma 
follows immediately from fl76|). Lemma is proved. 

It can be easily seen that, under the conditions formulated in Lemma 
4, for the functions n q (U), view ( [75] ) and (5C), also take place the relations 
similar to ( [751) , that is: 

Lemma 5. 

Under the formulated in Lemma 4 conditions: N = 2m + p, where p is the 
number of the annihilators of bracket (f|), on the submanifold A4, connected 
with the full family of m - phase solutions of (0), take place the relations: 

Corollary. 

On the submanifold A4' with the coordinates (J(X),9q(X)) take place the 
relations: 

{^(J(X)),^(J(F))} = eV^JAV], (82) 

{fc*(J(X)), n g (J(y))} = e 2 <[J, 6»S, e], - (83) 

where p a/3 ,v q are regular at e — > functionals of J(X) and 9q(X). 
Indeed, using ([74] ) and ([75]) , we obtain: 



{^(J(X)), fe^(JOO)} = e(E <(J(^))^(^)<5^ " n+ 

v=l OJ 

+ (*)) - Y)) + e*K[J, 91 e)^l(Y) = 

-e 2 K[J,9*,e]^l(Y). 

The relation flS3|) is proved by the same way. 

Now we shall calculate the brackets {J"(X), 9* a (Y)} and {9* a (X), 9^{Y)} 
on M.' . As was mentioned above, the functionals of type / q(X)J u (X)dX 
generate the local flows of form (pi]). The corresponding evolution of the 



26 



densities A u (ip, ip &a , . . .) of integrals a u (X), introduced in (|49| ) for the defini- 
tion of coordinates U U (X) in the vicinity of A4', also have the form similar 
to floTp, that is 



A%(9, X) = q(X)S u ^, (f ea ,..., e<f x , e<f 0aX , . . .) + 

+^QxS^{cp, 9e<*,..., etp x , etp gaX , ...) + ... . 
After the substitution of functions (f(8, X) in the form: 

<p(8, X) = ® in (0 + 6*(X) + - f X k(J(X'))dX', U[J, 9*, e]) (84) 

and integration with respect to 8, we obtain the expression for the evolution 
of a v (X) on the functions from Ai', characterized by the coordinates J (AT) 
and 8q(X), in the form: 

<»(x) = j:f; {k) [j,e*,q(x)](x)e k . 

k>0 

So that, the analogous expressions will be for the evolution of the func- 
tionals U(X) on the A4', since they are expressed in terms of sl(X) by the 
point substitutions. Besides that, since the system (^) generates at zero 
order of e on the functions from Ai f (having the form (|4])) the shift of the 
initial phases 8^ and does not touch the other parameters, the zero terms do 
not present in this series. So that, the Poisson brackets between U U (X) and 
J fi (Y) on the M.' will (as functions of coordinates J (A) and 9q(X)) have the 
form: 

{U"(X),J»(Y)}\ G{ e jX) = = Y,r; {k) [^e*]e k (85) 

k>l 

(that is {U V {X),J»{Y)}\ M , = 0(e) ate -> in the coordinates (J(A), #*(A)). 
Lemma 6. 

At Y 7^ X we have on the functions from A4', characterized by the coordi- 
nates J(X),8q(X), the relations: 

{TO/{K)} = a^[J l fi5,c] I (86) 

where a a/3 are regular at e — ► functionals of J and 8q. 
Proof. 

As can be easily seen, the functionals 8q ol {X) are expressed, according to their 
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definition (see ([56])), by the same (independent upon e) functions r a ($*(X)) 
in terms of the functionals t?* Q (X), introduced by the formula: 

x + - / A 'k(J(A"))«',UM(A'))<f"9, (87) 

e Jin 

by which the functionals 9q(X) are expressed in terms of $o(^0, introduced 
in flSTj). The formulas 

#*"(X) = r a (tf*(X)) (88) 

take place in the map on the A4', in which #o(X) takes at all e the same (that 
is independent upon e) values, that 6q(X) takes in the map corresponding to 
the functions r a in the definition of 9$ (X) in terms of t?o(X). 

The evolution of the functionals (^) according to the flows (|8"TD , gener- 
ated by the functionals of type / q(X) J^(X)dX, is determined by the evolu- 
tion of <f(9, X), given by the system (]61f) , and, besides that, by the evolution 
of J(X) and U[<^](X), which presents in $ in and in the coordinates J(X) 
and 9q(X) on the M! is given by the formulas ( |66"D and (85). After the sub- 



stitution of <f(9,X) in the form (84) into the system (|6lD and integration 
with respect to 9, removing the singularity at e — > in the shift of 9, we can, 
as in the case with U(X), conclude that the Poisson brackets of 9* a (X) with 



J^iY) on the A4' are regular (view (|88|)) at e — ► functionals of J(X) and 
9* (X). That is 

{9* a {X),J»(Y)} = '£e k \^[J,6Z). 

k>0 

Let now Y ^ Xq. Consider the functionals of type / q(Y')J fl (Y')dY' 
with the functions q{Y r ), having the support in the vicinity of Y, such that 
q(Xo) = 0. Taking into account the statement d67|) of Lemma 3, it can 
be easily seen, that the evolution of the functionals ( |3TD according to the 
corresponding flows has on the functions from the A4', characterized by the 
coordinates J and 9q, in the zero order of e at e — ► the following form: 

C(I)| £ .o = 7^- f 21 . • ■ / 2 '[?(I)Q^», ^(J)^, ■ • + 
{2n) m Jo Jo 

IX 



+®inj* (uP{J{X'))q{X')) xf dX> x ® m6ae0 ]d m 9 
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p2tt r27T 



r2ir r 2-n 

(after the substitution of tp{6,X) in form of ( |84D we use the relation ( |12"D for 
the functions QV)- The statement of the Lemma immediately follows from 
this. 

For the J(X ) (at the point X the depending upon e initial phase shift 
do is absent) we can obtain by the similar way on the Ai' (in the coordinates 
J(X),9* (X)): 

{9* a (X), J»(X )} = u«(J(X ))5(X - X ) + O(e), (89) 

(since at the zero order of e the functionals J q(X)J tI (X)dX generate the 
linear dependence of the initial phases 9q(X) upon the time with the frequen- 
cies q(X)u^(X) and at the point X : Oq(X ) = 9o(X )). For the functionals 
k a (J(X)), view (|75D, we have in the coordinates J(X),9%(X) on the M' at 
all Y, including X , the relations: 

{9* (X),k a (J(Y))}\ M , = O(e), e^O. (90) 



For the functionals ?9* a (X), introduced in fl87|), we obtain also, using ([82]) 
and (PD|), that their Poisson brackets with each other are regular on M.' in 
the coordinates J(X),9q(X) at e — > 0, and, so that, the same property is 
valid for the brackets of type: {9^ a (X), e o*( Y )}> that is 

{9* a (X),9t(Y)}\ M , = 1 a(3 [3,0*,e](X,Y), (91) 

where 7 a/3 (X, Y) are regular at e — > functional of J #q. We shall not need 
for the more precise information about the brackets of this kind. 

Let now consider the Dirac procedure of restriction of the bracket (fl5|) 
on the submanifold Ai', using the coordinates J(X),9q(X) and G l (9,X) in 
A 5 . 

For the Dirac restriction of the bracket (|45D on M' we must find for the 
functionals J U (X) and 6*Q a (X) the additions of type: 

V V {X) = . . jf* J v»(X, Y, 9, 6)0X9, Y)d m 9dY, (92) 



W a (X) = -i- fj. . . fj J wf(X, Y, 9, e)G\9, Y)d m 9dY, (93) 



(2tt)' 
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such that the flows, generated by the functionals J V {X) = J U (X) + V(X) 
and 0%{X) = Ql a {X) + W a (X), leave M' invariant. After that we must put 
on M'\ 

{r(x),j»(Y)y = {r(x)J»(Y)}\ M ,(j,e*), (94) 

{6* a (x),e?*(Y)y = mx)MY)}\ M >(J,e* ), (95) 
{r(x),e* a (Y)y = {r(x),e* a (Y)}\ M ,(j,9*). (96) 

On the functions v"(X, Y, 9, e) and wf (X, Y, 9, e) we obtain, respectively, 
the relations: 

x{G'(UM(I), <p{9, X), . . .), &(U[<p](Z), ip{ff, Z), . . .)}d m 9'dZ\ M , = 

= -{G ! (UM(I), p(9,X), . . .), rVp\{Y)}\ M < (97) 

x {G"(T%](X), <p(9, X), . . .), ff(UM)(Z), <p(#, Z),.. .)}d m 9'dZ\ M , = 

= -{G*(XJ[ V }(X), v (9,X),...),9* a [<p}(Y)}\ M/ (98) 

In the calculation of the Poisson bracket of constraints G l (9, X) with any 
functional on M! we can use the property, that on Ai' the values standing in 
the brackets in the definition of constraints G*(9, X) in terms of <£>(X) (see 
(|53"D) become zeros and, so that, in the calculation of the brackets of type 
{G l (9, X), C(9', Z)} on M.' we can omit the brackets of the kernels of the 
operators L, presenting in ([53"D, with C(9',Z) and replace the kernels of L 
from the Poisson brackets in the form of multipliers according to the Leibnitz 
identity. 

Regarding also the relations @, (g|), and (§§, flg§|), (0) for the 
functionals J and 9q, presenting in (|53[) , we can see that the Poisson brackets 
of type {G i (9 } X),G j (9',Z)} : {G\9,X), J U (Y)} and {G i (9, X),9* a (Y)} can 
on the submanifold M.' with the coordinates J(X),9q(X) be represented in 
the most general form at e — ► by the asymptotic series of type: 

{G l (9,X),G\9',Z)}\ M ,= 
= (2^F /o" ' 1^ £ Wll (' + o(X) + r + 9*(X) + ^1, x) x 
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X 



E ^ Up)i Jr + ^(X) + ^ x 

v n,fc>p>0 

s(Z) , , , s(Z) 



x ii[u[j,« v]] K + W + ^,C + W + ^^ ^^ m C, (99) 



s(X) and s(Z) denote the integrals k(J(X'))dX' and jf o k(J(X'))dX', 
and also: 

{G\9,X),r(Y)}\ M ,= 

= (2^p jf ■ • C z { e + w + ^r> T + ^ (x) + ^f 1 ' x ) x 

E ^ )>W (r + O *PO + [J, ^])^ W (^ " Y)) d m r, (100) 

^fc>p>0 e / 

{G*(0,x),r"(y)}U = 



C- ■ r Z W« ( d + W + r + O *(X) + x) x 



(27T) 

x f E r^ >w (r + es(x) + ^,[j,es])c**W(x-y)) <TV. (101) 

\fc>p>0 e / 

(All the functions in the sums within brackets, depending upon Z and Y, 
are replaced by the functions, depending upon X, according to the formulas 
of type: eS'(X - Z)f(Z) = ef(X)6'(X - Z) + ef x (X)6(X - Z). As can be 
easily seen, all differentiations with respect to X appear with the multiplier 

<•) 

The functions 

~«WKd e + *SP0 + 7 r k (JPO)^', X ) ( 102 ) 

C JA(] 

and 

*hmd 9 + ^ + 7 r k ( J ( x ')^', (los) 

are respectively left and right eigen vectors on the M.' of the linear op- 
erator in the space of 2tt— periodic with respect to 6 functions with the 
kernel {G l (6, X), G^{9\ Z)} and correspond to zero eigen values. The func- 
tions K[u[j t e* y€ ]](6 + 8q(X) + j Jx k(J(X'))<iX / , X), as can be easily seen from 
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( |100| ) and ( |101| ), are also orthogonal (at all e) to the right parts of (|9T|) , (|98|) , 
and, so that, the systems (l97[ ) , fl98|) are resolvable in the generic case. The 
solutions Vj(X, Y, 9) and w°(X, Y, 9) are defined modulo the arbitrary linear 

combination of vectors ft[u[J,6>*,e]](# + 9q(Y) + \ Ix k(J(X')dX' , Y), but it 
does not influence, view (|54l) , neither on the form of the additions V U (X) 
and W a (X), nor on the Dirac restriction of the bracket (f45|) on Ai' according 
to the formulas (p^)-(|9~6"l), since in the calculation of the brackets of V(X) 
and W a (X) with any functionals on the M 1 the contraction of this linear 
combination with the kernel of corresponding -£4[u[j,0*,e]](- • •) wm be replaced 
(according to the said above) from the Poisson bracket as a multiplier which 
is equal to zero on M.' . 

For the unique determination of the functions: 

v v 3 (X, Y, 9, e) = v»(X, Y,9 + 9* (Y) + - [ k(3(X'))dX', e) (104) 

and 

wf(X, Y, 9, e) = wf(X, Y,9 + 9*{Y) + - [ k(3(X'))dX', e) (105) 

in the coordinates J{X) and 9q(X) on M! we put N + m additional relations 
on them, demanding from them to be also orthogonal at all X and Y to the 
corresponding functions: 

+ 9*(Y) + - / y k(J(X'))dX',F,e), (106) 
e Jx 

that is 

■ ■ jf E v-iX, Y, 0, e)^ [Wo4] (9, Y, e)d™9 = 0, (107) 

j^i ■■■I E^(^^^^)^V[^oV]]^ y ' e )^ = ' ( 108 ) 

at all X, Y, e on M' . 

After the substitution of values v j(X, Y, 9, e) and u>"(X, Y, 9,e) in the 
form (1101,([105D into the systems (0), ©, where {(^(0, X), G»'(^, Z)}, 
{G i (9,X),J u (Y)} and X), 9* a (Y)} are taken in the form @, fllDTD 



K [U[J,0*,,]] 
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and ( |101| ) respectively, on the functions Vj and (after all differentiations 
with respect to X in (|97|) , (|98|) , appearing in every case with the multiplier e, 
the singular at e — > phase shift 9q(X) + i f x k(J(X'))dX', which presents 
in all functions, depending upon 9, can be omitted) we obtain the linear 
nonhomogeneous systems, which can be represented in form of regular at e — > 
asymptotic series with respect to e. As was mentioned above, in the generic 
case corresponding to the nonsingularity of the matrix of Poisson brackets of 
constraints in the finite-dimensional case, these systems are resolvable at all 
e, and, in the presence of the additional relations (|107|) , (|108|) , the functions 



v(X, Y, 6, e) and w(X, Y, 6, e) can be uniquely determined on the submanifold 
M' with the coordinates J(X), 0q(X) (the systems (p7[),(p8[) depend of them 
as of parameters) in the form of regular at e — > asymptotic series with 
respect to e, that is 

v](X, Y, 9, [J, 9*},e) = £ v v m (X, Y, 9, [J, 9* })e k , (109) 



fc>0 



\X, Y, 9, [J, 9*}, e) = £ w% k) (X, Y, 9, [J, 9* })e k . (110) 

fc>0 

Remark. 

Very frequently the Poisson brackets (|j) for the initial system (0), such as 
the Gardner - Zakharov - Faddeev bracket: 

{(p(x),ip(y)} = 8 r (x-y), 

or Magri - Lenard bracket: 

{<p{x), <p{y)} M L = -\ 8 "\ x -v) + (<P( X ) + f(y))S'{x - y), 

have such a form, that the corresponding to them bracket (45) is degenerated 
at zero order of e in the coordinates ip(9, X). This degeneracy arises because 
of the presence of the derivatives of 8— functions with respect to x (that is 
the operators d/dx) or the derivatives with respect to x of functions <p(x) 
in every term of such brackets, which arise in fl45|) with the multipliers e. 
However, in the introduced above coordinates J 1, {X),9q x {X) and G l (9,X) 
in the vicinity of M! the operators ed/ dX, being applied to the functions (ff5|) 
on the M, contain the nonvanishing at e — > value of type k a (J(X)) 
which permits to assume in the generic case (for the generic constraints) that 
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the systems on the functions Vj and Wj are nondegenerate in the zero order 
of e, and to write for the functions and wf the regular at e — > asymptotic 
expansions ( |109| ) and ( |110| ). 



More precisely, the expansions (|109|) and (|110|) have the form: 



v](X,Y,B,e)= £ v» {k)>(p) (X,e,[J,e*])e k 6M(X-Y), (111) 

fc>p>0 

w«(X,Y,6,e) = £ wf {kUp) (X,e,[3,e*})e k 5^(X~Y) (112) 

fc>p>0 

(here at every given k p runs the finite number of values p < k), the values 
V U (X) and W a (X), as can be easily seen, do not contain the generalized 
functions. The substitution of (|111|) and (|1 12|) in the corresponding systems 



(PTD,(P%D gives in the k-th order of e and at any given X and p (the coefficients 
with the derivatives of 5— functions must be equal to each other similar to the 
coefficients of powers of e) the linear (differential with respect to 9) nonhomo- 
geneous systems on the functions v v -^ ( p )(X, 0) and w®n.s ^,\(X, 9), depending 
of J and 8q as of the parameters. The right sides of these systems depend 
(in the linear manner) upon the previous V(k r ),(q) and W(k'),(q) k' < k. Under 
the assumption made above about the unique resolvability of the systems 
(|97D,(|98D at all e in the presence of the additional conditions ( pL07| ) and ( |108| ), 
corresponding to the nonsingularity of the matrix of Poisson brackets of con- 
straints on the Ai', the system described now will be uniquely resolvable 



under the additional conditions ( |107|) and ( |108|) in the corresponding order 



of e, and, so that, the series ( |111|) , (|112|) may be constructed by successive 
approximations. 

Besides that, as was mentioned above, the flows (0), generated by the 
functionals / q(X)J^(X)dX on the functions of form (0), leave invariant the 
submanifold M! at zero order of e, generating on it the linear dependence 
of the initial phases upon the times r' 4 , and, so that, in the expression ( |100| ) 
and, respectively, in the right side of the system (^) the zero term of e is 
absent [S^ = 0). From this we may conclude that, under the assumption of 
the unique resolvability of the systems on the functions Vj at zero order of 
e (in the coordinates J, 6*q on A4'): ^J(o),( p ) = 0, that is for the functions Vj 
and we may write: 

v»(X,Y,6,[J,6*},e) = 
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= E ^ {kUp) (x,e + e* (Y) + - f Y k(j(x'))dx f ,[j,e* ])e k d^(x-Y), 

k> P >i e Jx ° 

(113) 

wf(X,Y,0, M],<0 = 
= E + W + 7 / y k(J(a;'))rfX',[J,^])6^(X-r). 

fc>p>0 e Jx ° 

(114) 

The more precise information about the functions t>J and will be not 
necessary for us. 

Let us now formulate the main result of this paper. 

Theorem 1. 

Suppose that for the system Hamiltonian with respect to bracket (H) 
and having the family of m - phase solutions and the sufficient number of 
conservation laws take place all the described above properties of generic 
situation, concerning the functional independence of the parameters U u on 
this family and the possibility of the expression of parameters k, u and r in 
terms of them, and, besides that, takes place the relation, connecting the 
number of the annihilators of bracket (0) with the number of introduced 
above additional parameters r x ,...,r 9 (g = p). Then, under the assump- 
tion of regularity ((A) and (B)) of the introduced previously submanifolds 
M-Iqu k , (that is the possibility of constructing for any system (pD), generated 
by the functional / q(X)J u (X)dX , of the asymptotical solutions (|68| )), and 
the analogous regularity of A4 with the nonsingularity of matrix of Poisson 
brackets of constraints {G l (8, X), G^(6', Z)} on A4' at zero order of e in the 
coordinates J, 9q (that is the possibility of constructing of functions and 
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° in form of the asymptotical series ( |113| ),( |TT4"| )): 



1) The Dubrovin - Novikov bracket, defined by the formula (|3 



{IT(X), U"(Y)} = (A7)(U(X))5'(X -Y) + ^§^-(V(X))U x 6(X - Y), 



satisfies to the Jacobi identity. 

2) For the bracket (|38|) take place the following relations: 

{k a (U(X)), kP(U(Y))} = 0, {k a (U(X)), n q (V(Y))} = 0, (115) 
{k a (U(X)),U^Y))} =uZ(XJ(X))6'(X-Y)+^ x (V(X))6(X-Y). (116) 
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Proof. 

If the conditions, formulated above, are valid, we can restrict the Poisson 
bracket (15) on the submanifold A4' with the coordinates J(X),9q(X) ac- 



cording to the described procedure. According to the formula fl42|) , for the 
restricted on the M.' bracket we shall have the following formulas: 

{r(x),j»(Y)y = {r(x),r(Y)}\ M ,(j,9*)- 

r 2n r"2n 



(2< 



2m 







jj J j v»(X,X,9,e) x {G i (9,X),G^(9,Y)}x 
xvf{Y, Y, 9', e)dXdYd m 9d m 9', 



and the analogous formulas for {J U (X), 9* a (Y)}* and {6$*(X), 9f(Y)}* (with 
the replacing of functions v%, Vj to w^Wj for the functions 9^ a (X) and 
8q^(Y)). As previously, it can be easily seen that in the integration of 
functions (|TTB|) , QITi|) and |§) the singular at e -»• phase shift 0q(X) + 
i fx ~k(J(X'))dX', which presents in all functions, depending upon 9 and 9' 
(after the integration with respect to dX and dY), is unessential, and the 
Poisson bracket {...,...}* on M.' is regular at e — > in the coordinates J(X) 
and 9* (X). Using the relations (H, ©, ©, and O, (ED , we obtain 
for the bracket {...,...}* in the coordinates J, 6^: 

{r(x),j^Y)y = 

= e ({A7}(J(X))6'(X -Y) + (^(Q^)(J(X)))5(X - Y)) + 0(e 2 ), 

(117) 

{r(X),9* a (Y)}* = O(e) (118) 

at X ^ X , 

{05°(x) J ^(r)}* = o(i) (ii9) 

at e — > 0. 

^From this it can be seen that the Jacobi identities, containing only the 
functionals J at the points X, Y, Z, which do not coincide with Xo, coincide 
in the first nonvanishing order of e (at e 2 ) with the corresponding Jacobi 
identities for the bracket ( |3"8"D on the space of fields U U (X). Regarding now 
that the point Xq is arbitrary, and the expression (Rl) does not depend upon 
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Xq, we obtain that the bracket (|38| ) satisfies the Jacobi identity. Its skew- 
symmetry is a trivial corollary from the skew-symmetry of the bracket (H). 

The relations ( |1 15D ,( [TT5D now follow from (^), ( |83"D and ([7|) respectively. 

Theorem is proved. 

Theorem 2. 

Suppose now that under the assumptions formulated in Theorem 1 we have 
two different sets {I 1 , . . . , I N }, {I 1 , . . . , I N } of N conservation laws of system 
(0) of form (H), (some of the integrals of these two sets may coincide with 
each other). Then, the Dubrovin - Novikov brackets obtained with the aid of 
these two sets are coincide. That is, if U v = (V),U V = (P u ), where V \V V 
- are the densities of the integrals of the first and second sets respectively, 
(. . .) is the averaging on the family of m - phase solutions of (0), then the 
transformation of bracket (0), obtained with the aid of the set {J 1 , . . . , I N }, 
to the coordinates U V {X) on the M., expressed by the point substitutions: 

U" = tZ"(U) (120) 

in terms of U(X), gives the Dubrovin - Novikov bracket, obtained with the 
aid of the set {I 1 , . . . , 7^}. 
Proof. 

In the case of the generic situation under consideration, the Dirac restriction 
of the bracket (f45|) on the submanifold A4' is uniquely determined, that is, 
the brackets ( |117| ) - (|119[) , obtained in the coordinates (J, 8q) and (J, 9$) 



(corresponding to the first and second sets of integrals respectively) in the 
vicinity of A4', transform into each other under the corresponding transfor- 
mations of the coordinates (J, 9$) and (J, 6*q) on the M.' . Regarding (|5^) 
and the similar relation for tj and (J, #q), we can conclude that the transition 
from the coordinates J(X) to J{X) has the form: 

J V (X) = u u (3(X)) + £ k J(k)(JiJx, ■ ■ .,ikx,9l x , . . . ,0ofcx)> 

k>l 



and, so that, the transition from the brackets ( |117| ) to { J U (X), J tl (Y)} coin 



cides for such transformation (if X, Y X ) in the first nonvanishing order 
of e (at zero power of e) with the corresponding transformation of Dubrovin 
- Novikov bracket under the substitution ( |120| ), which proves the Theorem. 

The evolution of densities V v {x) according to the flows (|Tl"D , generated 
by the integrals 7 M , as can be easily seen from (|36|), has the form: 



V^, Vx ,...) = d x Q^( V ,^ x ,...). (121) 
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The flows, generated by the functionals / U^(X)dX on the space of fields 
U(X) with the aid of Dubrovin - Novikov bracket, has the form: 

U^ = d x (Q^)(XJ), (122) 

and, by such a way, represent (see (|30D ) the Whitham's equations for m - 
phase solutions of the Hamiltonian system, generated by the functional 1 11 (if 
it is not the momentum operator or annihilator of bracket (|3|) , let us remind 
that all these systems have the common family of m - phase solutions). 

All these flows commute with each other because of the commutation 
of functionals / J^(X)dX with respect to Dubrovin - Novikov bracket, and, 
besides that, the same is valid for the flows generated by the integrals with 
respect to X of the averaged densities of all functionals, having the form (||) 
and commuting with the Hamiltonian and integrals I u , since any of these 
functionals can be included into the set {I u } instead of any of presenting 
there integrals I u , and, according to Theorem 2, this will not change the 
bracket (|38|) . The integrals with respect to X of the averaged density of 
momentum operator (|5|) and the annihilators of bracket (^|), having the form 
(HD, generate in the bracket (|38D the shift with respect to X and zero flows 
respectively. The integral with respect to X of the averaged density of the 
Hamiltonian (^) generates Whitham's equations for m - phase solutions of 
system (0). 

In closing, the author expresses his gratitude to S.P.Novikov, who sug- 
gested the above problem, for his attention to the work and also to V.L.Alexeev, 
O.I.Mokhov, M.V.Pavlov and E.V.Ferapontov for fruitful discussions. 
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